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Abstract. We consider the effect of trapped branes on the evolution of a test brane 
' whose motion generates DBI inflation along a warped throat. The coupling between the 

■ inflationary brane and a trapped brane leads to the radiation of non-thermal particles 
on the trapped brane. We calculate the Gaussian spectrum of the radiated particles 
and their backreaction on the DBI motion of the inflationary brane. Radiation occurs 
either due to a parametric resonance when the interaction time is small compared 
to the Hubble time or a tachyonic resonance when the interaction time is large. In 

O |' both cases the motion of the inflationary brane after the interaction is governed by 

a chameleonic potential, which tends to slow it down. We find that a single trapped 
brane can hardly slow down a DBI inflaton whose fluctuations lead to the Cosmic 
Microwave Background spectrum. A more drastic effect is obtained when the DBI 
brane encounters a tightly spaced stack of trapped branes. 
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PACS numbers: 98.80.Cq, 98.70.Vc 
1. Introduction 

D-brane inflation [1, 2, 3, 4, 5, 6, 7] has become a time honoured subject whose 
development has led to plausible tests of string theory in a cosmological setting. With 
the launch of the Planck satellite [8] and the advent of precision cosmology, this 
possibility has become even more relevant [9, 10, 11]. Along these lines, there are 
three typical observables which may have an impact on models of inflation inspired or 
derived from string theory. The first one is the spectral index. The second one is the 
production of gravitational waves [12, 13]. Another one is the possibility of primordial 
non-Gaussianities [14]. It is quite likely that Planck will reduce the uncertainty on the 
spectral index in such a way that different scenarios such as brane inflation or modular 
inflation [15, 16, 17] may be distinguished. For instance, modular inflation models of 
the race track type tend to favour a rather low spectral index [18, 19, 20, 21]. Similarly, 
the absence of detectable gravity waves in these low field inflation models is generic. 
In this paper we will consider a related family of brane inflation models which do not 
fall within the category of slow roll models. In DBI inflation [22], a test brane evolves 
due to the steepness of the warp factor. This is particularly well motivated in string 
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compactification scenarios where a D3 brane evolves down a warped throat [23, 24] 
attracted towards an anti D3 brane. In general, potential terms such as a large mass 
prevent the existence of a slow roll period [6]. The brane may thus enter a regime of DBI 
inflation whose features are very different from the ones of slow roll inflation. In DBI 
inflation, the spectral index tends to be close to one [2, 14, 25] and non-gaussianities are 
large. Reheating in brane inflation [5, 26, 27, 28] is also quite different from reheating 
in standard inflation [29, 30, 31]. 

Recently, the issue of brane inflation in the D3 — D3 system has received a host of 
new developments. It has been realised that the potential does not only comprise the 
Coulomb interaction and a mass term responsible for the 77 problem [6] ; there are also 
corrections coming from bulk effects in the compactification scheme [32, 33, 34]. These 
bulk effects may require a detailed knowledge of the full compactification manifold. 
It happens that the leading contributions may have two effects. One possibility is the 
existence of an inflection point [35, 36] in the inflation potential when a fractional power 
of the inflaton in — c0 3 / 2 is present. Slow roll inflation is tunable in this context and 
requires a choice of initial conditions to overcome the overshooting problem [37, 38]. 
Another possibility is that the bulk term contributes as another mass term to the 
potential. In this second case, DBI inflation may occur. 

Our paper discusses brane models without attempting to realise their embedding 
in string theory. Hints along these lines can be found in [39]. In this paper, we will 
focus on DBI inflation as a field theoretical model corresponding to a warped throat 
with a quadratic potential. We will envisage the likely situation wherein the warped 
throat is not empty but contains a series of trapped branes [40]. These branes may 
be present due to fixed points in the compactification manifold. When the inflationary 
brane crosses these trapped branes, particles are created with a spectrum whose shape is 
Gaussian with a width related to the speed of the brane. The radiation of particles leads 
to the creation of a long range potential which might slow the motion down. The overall 
effect of the trapped brane is to induce a chameleonic term in the inflationary potential. 
Eventually this change of the potential implies that the Hubble rate is modified after 
the crossing of the stack. This effect decays rapidly as the radiated particles are diluted. 
As a result, the change of the Hubble rate is only effective on a few e- foldings after the 
interaction. 

In this paper we concentrate on the detailed analysis of the creation of particles on 
the trapped branes and the subsequent slowing down of the inflationary brane in the 
DBI context. In section 2, we recall some of the salient characteristics of DBI inflation. 
In section 3, we consider the creation of particles on the trapped brane. In section 4, we 
study the backreaction of radiated particules on the inflationary brane and the induced 
change in the potential. 
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2. DBI Inflation 

We are interested in brane models inspired from string theory. The main ingredients 
of our model will be the DBI nature of the inflaton dynamics and the quartic coupling 
of the inflaton to a scalar field \- The coupling is equivalent to the one between the 
inflaton and brane degrees of freedom on a trapped brane fixed along the inflationary 
valley [40]. In the slow roll regime, the same model has been considered in [39]. Here 
we analyse the case where the r\ problem prevents slow roll inflation but allows DBI 
inflation. 

The dynamics of the inflationary brane are described by the Dirac-Born-Infeld 
(DBI) action, where canonical kinetic terms have been replaced by a DBI term which 
can be expanded as a sum of higher order derivatives 

+V(<f>) + 9 ^x 2 \<P - 0i| 2 ) + / d'xV^^R (l) 

where R in the Einstein-Hilbert action is the Ricci scalar and is the metric. The 
metric is of the FRW type with no curvature 

ds 2 = -dt 2 + a 2 (t)dx 2 (2) 

in cosmic time. Conformal time is defined by a(t)di] = dt. We will exclusively focus 
on a warped case with an AdS throat where /(</>) = ^. The 't Hooft coupling 
A = R 4 /a' 2 = R 4 /lg depends on the radius of compactification R. In this scheme, 
the inflaton represents the radial distance between the moving D3 brane and the D3 
inside the throat : = VT^r where T 3 is the brane tension. The coupling constants g s 
and g are respectively the string coupling and the Yukawa coupling where g s ~ g 2 . 

The potential V(4>) consists of several terms. The Coulomb potential Vq3-D3 — 
D — 16 ^ 4 j with D = fd m ) i where / is evaluated at the tip of the throat, 
corresponds to the attraction between the D3 and D3 branes. There is also a mass 
term V2<f) 2 coming from radiative and supergravity corrections to the potential. In 
general the mass term V2 is positive apart from the case of a probe brane starting at 
the tip of the throat and moving towards the bulk [41, 42]. Corrections to the potential 
coming from bulk effects have to be added too. These corrections have integer and half- 
integer powers of and depend on the bulk of the compactification. The two leading 
corrections are proportional to 2 and 3 / 2 . In the first case, the model is inflationary 
when the branes are sufficiently apart under the influence of the quadratic potential. 
In the second case, the mass term and a negative 3//2 term imply the existence of an 
inflection point around which the potential becomes V + Vi(<f> — inflection) ■ 

There is also a trapped brane whose location is fixed at 0i along the throat. 
Particles \ 011 t ne trapped brane are coupled to the inflationary brane with a quartic 
coupling at leading order. We will see that this coupling is responsible for the slowing 
down of the inflationary brane when crossing the trapped brane. 
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We define the equivalent of the Lorentz factor 7 

7 = y/i + f(<p)g, u d»<m = ^T^r (3) 

The inflationary dynamics are governed by the Klein Gordon equation for the spatially 
homogeneous inflaton field 

l^wX^X+ ldV 1 # 1 d f n (A\ 



3H x , 1 , 1 d/ \, 1 d/ , 3 i2 l d/ 



or equivalently 

^ + ^ 0+ ^l : d0 + ^J-^ + i 02 75 = o (5) 

In the ultrarelativistic limit, both the potential term and the friction term coming from 
the expansion of the universe are subdominant. 

With the warped factor /(0) = ^, equation (4) becomes : 

A,M^M^^M ldV # 3 40 3 20 2 

+ + -0 +-—- — + — + — =0 6 

7 7 d(p A7 A7 2 

It is easier to analyse the DBI dynamics using the Hamilton- Jacobi formalism. This 
formalism consists in eliminating the dependence in t and using instead 0-dependent 
functions. It is valid provided <f>(t) does not oscillate. Using the Friedmann equation 

3MJ, 3M p9s V / 
where Mp is the Planck mass, the potential can be written as 



V(4.) = 3M 2 p g,H 2 -7 + 7 (8) 

If we differentiate (7) and upon using the Klein-Gordon equation we can express ^ and 
obtain 

<P=-2M 2 P g s -— (9) 

7 a<f> 

Similarly, the 7 factor can be written in terms of functions of 

7(0) = ^l + 4M^/(0) J^j (10) 

For a given potential, the Hubble function of is the solution of a differential equation. 
Focusing on the limiting case where — > as £ — > 00 and expanding H(<j>) in powers of 
0, the relativistic factor is dominated by 

dH (ii) 



7(0) « 2M%g,*JW) 



d0 

As a result, the inflaton follows a universal trajectory dictated by the coupling /(0) 

^"TM a,ld ^ (t) = x (12) 
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The potential term has an influence on the scale factor only. 
We will exclusively focus on a quadratic potential where 

V{4>) = m 2 <p 2 (13) 

In this case, there are co-dominant terms in (8) and one finds 




3v/A \ V M p9s 
The scale factor becomes 
a I t x 



where 



a V £ 




A m 



(15) 



(16) 



3g s M P 

Inflation occurs when e < 1. We will focus on the case where e « 1 which leads to 
interesting restrictions on the parameter space (see section 4). In all cases the coupling 
— and jj^ can be adjusted to verify e < 1 at strong coupling A >> 1. The 7 factor 
becomes 

7 = (it) 

A e 

In the following we shall use conformal time to study the creation of particles when the 
DBI brane crosses the trapped brane. In conformal time we find that 



in such a way that 

H « -?T\ — ~ 2r? _2 (19) 



a" 



We deduce 6 oc (— ) m (^Y to leading order. We will make use of these identities 
in conformal time when discussing particle creation. 



3. Particle Creation 

3.1. WKB approximation 

So far we have not taken into account the presence of the trapped brane. In fact, the 
trapped brane has an influence on the inflationary brane evolution. To determine its 
effect we need to study the quantum modes of the field \. Let us expand the quantum 
field x i n terms of creation and annihilation operators. For convenience sake we work 
in conformal time. The field reads 

/c\ 3 k 
j^ 3 (a kXk (v)e ikx + alxl(v)e- ikx ) (20) 
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and each mode satisfies the Klein-Gordon equation 

Xl + 2H X ' k + k 2 Xk + a 2 g 2 \<P - fa\ 2 x k = 



We can put in a Schrodinger form by defining 

v(v) = a x(v) 

So \1/ obeys 

with a time dependent frequency 

Mv) = Vk 2 + A( V ), A( V )=a 2 g 2 \ 

where is the unperturbed 4>(rj) corresponding to the unperturbed motion of the 
inflationary brane. Here we have assumed uj 2 positive. When it is negative, the regime 
is tachyonic and the frequency will be written Vt k = ±ic<Jk = ±i-\/|k 2 + A(?7)|, so that 
n 2 k = -co 2 > 0. 

The equation for the modes (23) can be approximately solved using the Wentzel- 
Kramers-Brillouin or WKB approximation. Far in the past, the solution is assumed to 
be in a Bunch-Davies vacuum where for 77 — > —00 we have 



a 



(21) 
(22) 
(23) 

(24) 



1 



y/2u k (ri) 

When the inflationary brane has passed through the trapped brane and rj 
mixture of two possible modes 



(25) 
0, it is a 

(26) 



y/2uk(v) \/ 2u; k(r}) 
where a k and /3k are the Bogoliubov coefficients. The WKB approximation is valid when 
1 41 < 1- We define 



and find 



R 



R = 



UJ 



n g 2 a 2 \4>-<t>i 


2 + <j)'g 2 a 2 \ 


0-01 


- \ (f )' 


(P + aV|0-0i 


2 _ sly/2 



(27) 



(28) 



Our goal is to determine under which conditions the WKB approximation is 
violated. There are two physically different situations where the analysis can be easily 
carried out. They depend on 

- " 2 (29) 



9W 

which is a constant in DBI inflation with a quadratic potential 



1 



m 



ge 2 V\ 3g 3 \M Pi 
The creation of particles is very different for small or large £. 



(30) 
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3.2. Small £ behaviour 

In this region of the parameter space, the creation of particles occurs when the DBI 
brane is close to the trapped brane which corresponds to the regime where 

(31) 

In this case the R factor becomes 

|0Va 2 |0-0!| - \ , 
~ ncar ~ |F + aV|0 - 0i| 2 " vl 3/2 
and non-adiabaticity is mainly present in the region where g 2 a 2 \(f) — 0i| 2 ^> ^ (or 
equivalently g 2 |0 — 0i| 2 3> 2i/~ 2 ) for which we have 



R 





0i 1 


(F + aV|0- 


^| 2)3/2 



(33) 



From now on we define 



/C = - (34) 
a 

which corresponds to physical momenta, 

The creation of particles arises when R > 1 at its maximum which is located at 
g 2 \4> — 0i | 2 = ^ and is given by 

R ~ — ^ (35) 



implying that the creation of particles occurs when /C < ^ry5'|0| • The maximal 
extension of the non-adiabatic region is 

(36) 




9 

Notice that A0 > - as £ < 1. 

Inside a region of size |0 — 0i| < around the origin, there is a domain of 

non-adiabaticity when 

2\k 2 - ^| 3 / 2 



l(-)'l 

^ nl ,°.L, >l (37) 



a 



corresponding to 

V^2 - 2 2 / 3 # < /C < x/2 + 2 2 / 3 # (38) 

Due to the very small width of this zone in momentum space, virtually no particles are 
created in this interval. 

Finally when /C < V2H, there is a tachyonic instability. The size of the tachyonic 
region is much smaller than the size of the non-adiabaticity region where most particles 
are created. In fact, the creation of particles in the tachyonic region is negligible as it 
scales like exp(if At tachyon ) ~ exp(£) ~ 1 where the time spent in the tachyonic region 
is At tachyon ~ ^tp 
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It is interesting to notice that the time spent by the brane in the interaction region 



is 



At 



H 



(39) 



implying that the interaction lasts less than a Hubble time and therefore the interaction 
time corresponds to a number of e-folds 



AN = — « < 1 

a 



(40) 



The interaction is almost instantaneous. 

So far we have used the fact that the Hubble rate is nearly constant in the interaction 
region. This can be ascertained as the variation of the Hubble rate in the interaction 
region is given by 



AH 



H 



1 



(41) 

Therefore we must impose that gV\ 3> 1. We shall see that this is always the case 
when e<l. 



3.3. Large £ behaviour 

In this case, the ratio R can be simplified in two regimes depending on whether the DBI 
brane has moved far from the trapped brane or not. If the inflationary brane is far from 
the trapped brane then \<f> — <pi\ ^> <p/H and 



|^ 2 a 2 |0-0x 


2 _ 


1 (a" 

2 \ a 


)'l 


\k 2 + a 2 g 2 \(p — 


01 


2 _aT_ 

a 


3/2 



R « i? far = — — 2 v % ' (42) 

|A; 2 + aVl0-0i| 2 - vl 3/2 

On the contrary, if they are close, \<p — 4>\\ <C <j)/H so 

R~R \iSl\ (A3) 

ncar |F + aVl0-0i| 2 - vl 3/2 ~ 2 l fc2 - vl 3/2 
These two regimes capture all the physics of the x-particle creation. 

First, we will study the case when the two branes are far from each other. The 
ratio R is maximal when g 2 \4> — <pi\ 2 = 2/C 2 + (2 — 7e)H 2 and its value is simply 



\2HIC 2 


-4eH 3 \ 


|3/C 2 — 


6e #2|3/2 



-Rmax — , Qjr2 CTT2I3/2 



We find that there is a pole for /C po i c = \2eH . The value at the origin is very large: 
-R m ax(^C = 0) k e -1 / 2 > 1 as e < 1. In all the interval between the origin and the 
pole, -R max is greater than one. We want to determine the physical momentum /C max for 
which i? max = 1 and then becomes smaller than unity for larger momenta. Expanding 
around the pole 

2 H 

^ " 37K^-^ ol e) 1/2 { } 
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implying that 

AH 2 

Cax « -27T- (46) 

Hence we find that for physical momenta < /C < /C max ~ the WKB approximation 
is violated around g 2 \<p — <pi\ 2 = 2(/C 2 + H 2 ). Therefore there are two non-adiabatic 
regions far from the trapped brane: regions I and II centered respectively around <pA and 
(f>B- The approximation used here is valid as the maximal extension of the interaction 
zone is given by A0 y implying that A0 ^> when £ 1. 

We now consider when the inflationary brane and the trapped brane are close to 
each other. From (43), we find that R > 1 when 



'2 - 2 2 / 3 H < JC < V 2 + 22/3/7 (47) 

The inequality (47) gives the range of physical momenta for which a non-adiabatic 
region appears in the immediate vicinity of the trapped brane. It is consistent with 
# 2 |0-0i| 2 <2# 2 as£>l. 

On top of the non-adiabatic instability detailed above, there is a tachyonic resonance 
when uj\ < corresponding to 

K? + g*\<i> - fa\ 2 - 2H 2 < (48) 

So for physical momenta larger than 2H 2 , there is no tachyonic regime. We define 7]_ 
and i] + as the two turning points such that = ujI(i] + ) = 0. There is a physical 

momentum JC* for which the non-adiabatic and tachyonic regions intersect in just one 
point. For < /C < /C*, the regions intersect and for JC* < /C < /C max the tachyonic 
region and the non-adiabatic region are disconnected. In the tachyonic regime, we also 



use the WKB approximation. It is valid when 



< 1 and the modes are 



The Bogoliubov coefficients change after going through a non-adiabatic region. 
The interaction time can be estimated and gives 

HAt w ^ (50) 

leading to 

AN = — w e (51) 

So the interaction region is spread out over a large number of efoldings. Moreover we 
must impose that the variation of the Hubble rate is small in the interaction region 

AH ~ 1 (52) 

-ir-^Tx (52) 

hence we must have gey/X ^> 1. 
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3-4- Creation of particles 

We are interested in the particles created when the inflationary brane crosses the 
trapped brane. This happens when the WKB approximation breaks down. Let us 
first concentrate on the £ <C 1 regime. The calculation of the number of particles in this 
region is well-known [40]. The Bogoliubov coefficient is obtained by expanding 



around <p\ (5t = t — t\) and integrating around a contour where the WKB approximation 
is still valid in the complex plane. The end result is that 



where /Ci = kfa\. As a result, the spectrum is Gaussian with a width determined by 
the speed of the brane. 

In the £ 3> 1 regime we have found different situations depending on the physical 
momentum /C. The configuration where < /C < /C max is the most complex. There 
are two symmetric non-adiabatic regions far from <pi and a tachyonic region in the 
proximity of <\>\. The tachyonic region intersect with the non-adiabatic regions for 
< /C < fC* < JC max . For /C max < /C < \Jl — 2 2 / 3 if, the tachyonic region is still there 
but we no longer have any non-adiabatic region. For v2 — 2 2 / 3 H < /C < V2H, there is 
a tachyonic region and inside of it a non-adiabatic zone. For y/2H <JC< V2 + 2 2 / 3 H, 
there is no tachyonic resonance but there is a non-adiabatic region around <j)\. And 
finally for any K, > \Jl + 2 2 / 3 H, there is no tachyonic resonance and the regime is 
always adiabatic. Notice that all these intervals are easy to interpret using the physical 
wave number K. The size of the physical intervals in K, is time- independent. 

Let us first study the creation of particles for < fC* < fC < /C max - This case is 
typical and will allow us to deduce the particle spectrum in the other intervals too. 




(53) 




(54) 




tachyonic zone 




— OO 



non-adiabatic I 




non-adiabatic II 



Figure 1. Configuration of the interaction zone in the complex plane 
for K,* < K. < /C max 
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Initially, for rj — > — oo, the modes are in a Bunch-Davies vacuum. Then the WKB 
approximation breaks down in the non-adiabatic region I, far enough from it the modes 
become 

^(q) = c-ifuMW , ^(^) e iPu, k ( v ')d v ' ^ 

We need to determine the two Bogoliubov coefficients. In fact they can be expressed as 
transmission and reflection coefficients i and r 

«fc(^ < Va) = ta k (r) > 7] A ) and /3 k (r) > r) A ) = ra k (r] > r) A ) (56) 

where the conformal time 77 increases. So with the chosen initial condition a k (— 00) = 1, 
this gives 

ot k (v > Va) = l/t, (3 k (v > Va) = r/t (57) 

The solution before the non-adiabatic region I is linked to the solution on the other 
side of the non-adiabatic region by an analytic continuation in the complex plane. We 
consider that the time variable 77 is a complex variable. We draw a contour around 
7] — t]a in the complex plane (see figure 1). The radius of this semi-circle must be 
large enough for the WKB approximation to be valid before and after the non-adiabatic 
region. This is the case if 

Wk(VA)\V ~ Va\ > R(Va) = ^max (58) 

and we trust the expansion of u^ 1 around tja 

1 1 ful\ , . Afi M 

Vk(v) v k (r} A ) 

implies that 



-f (VA)(V ~Va)-J2 -?T-^ - ^ (59) 
where is the n-th derivative of R at i]a- The condition on the size of the semi-circle 



Mv) ~ - p — r V ( 60 ) 

RmzxiV ~ Va) 

Using the fact that in the non-adiabatic region uj k (r] A ) = O(aH) and Ra = O(^) = 
0(1), the condition on the contour is equivalent to 5t > 1/H. Now the duration of the 
interaction region is implying that the contour circling around the non-adiabatic 
region is much smaller than the interaction region. As a result we obtain that 

exp^/^w) «-p( T i/ fl _ (r J,_ r)4) d,/) (61) 

Notice that 77 — 77^4 is negative before the non-adiabatic region and positive later on. A 
positive-frequency mode is changed in a negative-frequency mode when going through 
the non-adiabatic region (and vice versa). With (77 — r] A ) = pe ld and di] = ipe l6 d9, 
the contour integral is given by the residue theorem. A factor of i also appears from 
Tj-rjA^ e m (r] A - v) in 7^= so finally 



(3 k = -iexp(-2i£ A ) exp (jL-^j = -iexp(-2i£ A ) exp {^^^ 



(62) 
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with Ka = k/a A and the phase is 

/VA 
co k drj (63) 
-oo 

Probability conservation imposes |o;fc| 2 — \fik\ 2 = 1, so we deduce 



a k = e lip ■ 



'l + exp [sVZttj^) (64) 



where ip is a random phase. 

We now compute the Bogoliubov coefficient b k of the non-vanishing wave in the 
tachyonic region [43]. We neglect the decaying mode in this region depending on the 
coefficient a k . We draw a contour around the turning point r]_ where u\(ri_) = 0. We 
assume the radius of the semi-circle is large enough for the WKB aproximation to be 
valid along the contour : 

"1 = ^(V-)(V-V-) (65) 

After an analytic continuation we find 

fe = a»e^ + »+&eK*- + *> (66) 
with the phase 

9~= f ujkdr] (67) 

J — oo 

depending on the wave evolution before the turning point. Then passing around the 
second turning point at we obtain a new contribution to the Bogoliubov coefficients. 

h = e J " nd V , ('- + *>6 fc = e #- ndv (ft + a k e~^-^ >) (68) 

a fc = ^-+i) 6fc = e'- ^ (a fc + fte 2 ^*)) (69) 

The final contribution comes from the non-adiabatic region II where the wave appears 
to break into reflected and transmitted ones 

a kX ~ ^R- X+ + T^x- (70) 

PUX+-+R+X- + T+X+ (71) 
We are particularly interested in the the final Bogoliubov coefficient 

Pi = R.+T + (72) 
The first situation (70) is the same as in (56) but with a phase 9 B . 

R- = T -a k = -iexp(-2i# B ) exp (^j 1 ^ j «k (73) 

where Kb = k/as The second situation (71) is the dual configuration for where we 
find 

* 2 



\T + \ 2 = \p k \ 2 1 



(74) 
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via current conservation. Therefore we can introduce another phase d such that 



T + = ^fa\]l + e^^ (75) 
As a result we find the Bogoliubov coefficient 

Let us now examine the other regions. For /C 2 nax < /C 2 < (2 — 2 2 ^ 3 )H 2 , there is a 
purely tachyonic contribution. For a.k{— oo) = 1, we know from (68) that 

|^| 2 = e 2/ -^ (77) 

The same behaviour is valid for (2 — 2 2 ^ 3 )H 2 < K? < 2H 2 because the non-adiabatic 
region has no importance when it is included in the tachyonic regime. Indeed the passage 
through the non-adiabatic region in a tachyonic region only changes the phase of the 
Bogoliubov coefficients. 

When 2H 2 < K? < (2 + 2 2 ^ 3 )H 2 , there is a non-adiabatic region centered at (f) 1 . 
In this region where /C is large, the duration of the interaction is less than 1/H. We 
continue the wave function in the complex plane on a semi-circle of radius larger than 
1/H such that 

— « J !C 2 - 2H 2 + g 2 <p 2 5t 2 w gfaSt + ^ -. ^ (78) 
a v 2g<p 1 5t 

Notice that the WKB approximation is valid along this circle. The integral J ^dt picks 
up an imaginary value due to the residue at the origin. As in (62) positive and negative- 
frequency modes are exchanged. And the coefficient (3^ is computed with the residue 
theorem 

2H 2 -K 1 2 

\f3{\ 2 = e (79) 

Notice that this is a small number in the region where K > 0{y/2H). 

For < /C < /C*, the method used for JC* < /C < /C max is not reproducible as the 
semi-circles drawn around the non-adiabatic regions I and II intersect the semi-circles 
drawn around the turning points. We will see that the low energy part of the spectrum 
is dominated by the tachyonic instability and therefore does not really depend on the 
non-adiabatic region. As a result we will extend our results to the whole momentum 
range < K, < /C max . 

We have now obtained the Bogoliubov coefficients in both regimes for £. This is 
particularly important as one can define an adiabatic invariant for equation (23) 

Nk = ^(^ + \*A-\ = \^ (80) 



2 \ u 2 ' J 2 

As explained in [30], this is the comoving occupation number of particles with 
momentum k. 
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4. A DBI Chameleon 

4-1. The modified potential 

We will compute the energy density of the created particles. In fact this energy density 
appears to add a contribution to the potential of the inflaton. This effect corresponds 
to the slowing down of the inflationary brane by the trapped brane. The emitted energy 
density is 

Px = J ^ UkUk (81) 

where the integration is done over the physical momenta. This brings an extra factor 
of 1/a 3 corresponding to the dilution of the created particles. The frequency Cj^ = ui^/a 
is the rescaled frequency coming from the energy in conformal time |cj^|x| 2 = 
Far from the trapped brane, once particles have been created in the immediate vicinity 
of the brane, 

ss— yi&i (2^)s (82) 

This approximation is valid as long as g\<j> — <pi\ ^> y^\2H 2 — K, 2 \. When ^ < 1, the 

right hand side is at most \J 2g\(j)\, which implies that |0 — 0i| 3> i.e. outside of 
the interaction region. In the £ >> 1, the same condition gives \<f> — 0i| 3> y which is 
also the size of the interaction region. In the effective description of the inflationary 
brane motion, this energy density is equivalent to a linear potential and a constant force 
towards the trapped brane. In the absence of the inflationary potential, this force may 
pull the passing brane towards the trapped brane. We are in a position to determine 
the effective potential in both cases £ 1 and £ ^> 1. 

When ( < 1 the potential can be easily calculated as the particle spectrum is 
gaussian 

p^7^%(g\4>i\f /2 9\4>-^\ (83) 

where as = ai and the energy density is diluted before finally tending to zero rapidly. 
In the £ 3> 1 regime, the creation of particles is largely dominated by the tachyonic 
instability in the vicinity of the trapped brane and depends on f* ^fdt. The behaviour 
of this integral is dominated by the region of small momenta. Let us concentrate on the 
low frequency part of the spectrum. In this case the integration region is maximal and 
the turning points are located at 

8t± w ±V2^- (84) 

#l0l| 

whose norm is always less than t\. When K is small, we can expand VLk 

^^^2H 2 -g^\5V- 1 =- 2 (85) 

a 2^/2H 2 - g 2 4>\8t 2 a 
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Integrating between the two turning points we find 

t+ tt k , 7t 2H 2 -Kl 

—dt « = (86) 

_ a 2 

where /Cg = k/as and we have defined 

dw 1 



1 1 1 r 1 



v / r^?(i + ^ M )2A 



(87) 



where we recall that ge\f\ ^> 1. This result is very important as it shows that the 



spectrum is approximately Gaussian with a width y g\(fii\ which is much smaller than 
the band of integration A/C = y/2H. Moreover the amplitude of the number of particles 
depends on e 2nH2 ^ 9 ^^ which is an exponentially large number. Notice that the same 
factor is negligible in the ( C 1 regime. Now that we know that the spectrum is 
dominated by the tachyonic instability, we can evaluate the energy density 

^-T^^ 2 ^ 11 ^!^!) 37 ^^-^! (88) 

Notice that the main difference with the £ <C 1 is the presence of a large exponential 
factor coming from the tachyonic instability close to the trapped brane. Here also, the 
energy is diluted after going out of the interaction region. The scale factor a s takes into 
account the fact that momenta are red-shifted in the integration region 

In both cases the effective potential after the interaction region becomes 

V ~ m 2 2 + -J^y^H 3 ^ - 0x| (89) 

where the coupling function depends on £ and reads 

1/(0 « r 3/2 (90) 

when £ <^ 1 and 

V(0 « 9r 3/ V^ (91) 

when £ ^> 1. This potential is similar to the ones used in chameleon models [44]. It has 
a moving minimum where 

This minimum goes to the origin as the scale factor increases. This implies that the 
effect of the trapped brane is only relevant for a few e-foldings after the passage through 
the interaction region. Immediately after the passage, the minimum is located at 

0! 16vr 3 3# Mj v ; 

This gives a criterion for the influence of the minimum of the trapped brane on the 
motion of the inflationary brane. If 0j n i <f>i, the trapped brane has no influence as the 
inflaton feels the m 2 2 branch of the potential. On the other hand, if H is large enough 



Brane Bremsstrahlung in DBI Inflation 



16 



and 0i n i > 0i, the inflationary brane feels the steep potential due to the trapped brane. 
In this case, the motion of the inflationary brane is affected for a few e-foldings. 
Let us evaluate the jump in the potential at the end of the interaction zone 

^ ~ 7^r 3/2 e 2ni H 4 (94) 
(27r) 3 

when £ 1 and 

AV " (2W < 95 > 

when ^Cl. In both case it depends only on H 4 and £. For large enough H, this jump 
in V can also change the Hubble rate due to the release of energy in the form of radiated 
particles. 

As we have considered that H is nearly constant and we have neglected any backreaction 
on the dynamics of the inflaton while the brane particles are created, we must impose 
that AV/V < 1, which gives for £ < 1 

^ < 3. WW (96) 

p 

And for £ > 1, 

^ < 3.( 27 r)W/V^ (97) 
p 

This condition gives an upper limit on £. In fact, if we require that the Hubble rate H 
should be at least of order 1 GeV, £ must be at most equal to £n m rs 14 for g ~ 1CT 1 . For 
larger values of £, the tachyonic instability implies that there is a strong backreaction 
due to the explosive creation of particules. In this case, the inflaton loses all its energy 
very quickly and transfers it into radiated particles. 

4-2. Discussion 

Let us now make explicit some of the constraints on the parameters of the model. Using 
the condition e < 1, we have 

m > -^=M p (98) 
v A 

Hence the masses must be large enough and related to the Planck mass. This implies 
that 

{ y> irx (99) 

guaranteeing the constancy of H in the £ 1 region. Using this lower bound we obtain 
that 

a 2 

m > ^M p (100) 

For a fixed string coupling g s , we see that £ determines the range of masses leading to 
inflation. 
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and 



When £ is large, the mass m can be smaller than the Planck scale. In the large £ 
regime we have 

i » w p < 101 > 

A > — (102) 
m 

This regime tends to favour small masses and a large compactification radius. 

So far we have not tried to connect DBI inflation to observations. Let us now 
assume that the DBI inflation regime we have analysed is responsible for the phase of 
inflation resulting in the quantum fluctuations leading to the CMB spectrum . In this 
case, the COBE normalisation determines the curvature perturbation 

So £ ~ g 2 £ ~ 1CT 5 . For a reasonable value of the string coupling g s ~ 1CT 2 , we find that 
inflationary branes whose quantum fluctuations lead to the CMB anisotropies must be in 
the £ 1 regime. We have another observational constraint ; since gravitational waves 
have not been detected yet, the ratio r — ^ must be small. The tensor perturbations 
spectrum being 

n - 1| d04) 



We deduce 



16e 

r = — < 1 (105) 

7 

With both (103) and (105), we find 

— < 10~ 5 (106) 
M p ~ 

We note that this upper bound for the Hubble rate is smaller than the bound given by 

(96) and therefore the backreaction is always negligible in the £ <^ 1 regime. 

We have seen that in this regime the slowing down by a trapped brane can be effective 

5^ 



or equivalently 



l| > * (108, 



Using (100), we find 



^ > 16v/3VV 1/2 (109) 

p 
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This is not compatible with (106). So it is not possible to slow the motion of the 
inflationary brane drastically after hitting a single trapped brane. If we assume that 
there exists a stack of N closely packed branes in the interaction region, then equation 
(93) becomes 

„ 1 iHiWA^ * 
0! 16tt 3 3g Ml y ! 

And condition (109) becomes : 

llGv^C 172 (HI) 



#2 



Unless we have at least N ~ 10 9 branes in the stack, the slowing effects of the stack is 
not drastic. The motion of an inflationary brane leading to the CMB spectrum is hardly 
affected by trapped branes. On the contrary, branes in a regime £ 3> 1 are efficiently 
slowed down because of radiated particles: a brane bremsstrahlung. 



5. Conclusion 



We have studied the slowing down of an inflationary brane in the DBI regime along an 
AdS throat when it hits a trapped brane. We have shown that the brane motion may 
be slowed down by the creation of particles on the trapped brane. This creation occurs 
either by parametric resonance or tachyonic instability. The latter case happens when 
the interaction region is large compared to the Hubble rate. Once the brane has left the 
interaction region, the effect of the radiated particles is to generate a linear potential 
whose slope is greatly enhanced in the tachyonic case. This is enough to stop the brane 
for a few e-foldings until the number of created particles has been diluted. Branes 
crossing the trapped brane in a time smaller than a Hubble time, are not drastically 
affected by the radiation of particles as they are slowed down in manner which does 
not alter their motion very much. A dramatic effect is only possible when a very large 
number of trapped branes are stuck in the interaction region. 

An interesting consequence of our results is a selection mechanism for the motion 
of branes in a throat where trapped branes are also present. Indeed, as the branes 
go past the trapped branes, the ones leading to a tachyonic instability are successively 
slowed down. Only the ones with little interaction with the trapped branes can keep 
on moving. It happens that the ones with small curvature perturbations are the ones 
which would be little hampered by the presence of trapped branes. This could have 
an interesting effect selecting the motion of branes whose inflaton would generate small 
levels of curvature perturbations. 

Another interesting possibilities is the modification of the perturbation spectra as 
the DBI brane crosses one or many trapped branes. Even if the motion of the DBI brane 
is not altered, the slowing down term in the potential having an effective relevance during 
a short time, one may envisage that the spectrum of primordial fluctuations might be 
affected by a jump in the derivative of the potential. This may even be the case when 
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the slowing down of the inflaton is negligible. A thorough study of this possibility as 
well as models with a jump in the potential energy is left for future work. 
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